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ABSTRACT 



What are the frequency, shape, kinematics, and luminosity of Lya envelopes surrounding radio-quiet quasars at high redshift, and 
is the luminosity of these envelopes related to that of the quasar or not? As a first step towards answering these questions, we 
have searched for Lyce envelopes around six radio-quiet quasars at z ~ 4.5, using deep spectra taken with the FORS2 spectrograph 
attached to the UT1 of the Very Large Telescope (VLT). Using the multi-slit mode allows us to observe several point spread function 
stars simultaneously with the quasar, and to remove the point-like emission from the quasar, unveiling the faint underlying Lya 
envelope with unprecedented depth. An envelope is detected around four of the six quasars, which suggests that these envelopes 
are very frequent. Their diameter varies in the range 26 < d < 64 kpc, their surface brightness in the range 3 x 10~ 19 < p. < 
2 x 10~ 17 ergs -1 cnr 2 arcsec~ 2 , and their luminosity in the range 10 42 < L(Lya) < 10 44 ergs -1 . Their shape may be strongly 
asymmetric. The Lya emission line full width at half maximum (FWHM) is 900 < FWHM < 2200 kms~'and its luminosity 
correlates with that of the broad line region (BLR) of the quasar, with the notable exception of BR2237-0607, the brightest object in 
our sample. The same holds for the relation between the envelope Lyo- luminosity and the ionizing luminosity of the quasar. While 
the deep slit spectroscopy presented in this paper is very efficient at detecting very faint Lyo- envelopes, narrow-band imaging is 
now needed to measure accurately their spatial extent, radial luminosity profile, and total luminosity. These observables are crucial 
to help us discriminate between the three possible radiation processes responsible for the envelope emission: (i) cold accretion, (ii) 
fluorescence induced by the quasar, and (iii) scattering of the BLR photons by cool gas. 
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1. Introduction 



galaxies: quasars: emission lines - galaxies: quasars: individual: SDSS J14472+0401, 



The so-called "Lya blobs" have attracted much attention in 
the past few years. They are extended nebulae at high redshift 
(z ~ 2-5) emitting in the Lya line, with typical sizes of 
~ 10" or ~ 100 kpc Jsteidelet dJ 12005 iMatsuda etal]|2004 
1201 It I Yang etail 120091 7 Their Lya luminosity is typically ~ 
10 43 erg s , and, while they are numerous at hi gh redshift, they 
seem to disappear at moderate to low redshift dKeel et al.l l2009). 
The source of their power has remained unclear and a contro- 
versial subject. Some host an active galactic nucleus (AGN) at 
their centre, whi le others apparently do not dGeach et al.| [2009; 
Yang et al. 2009). Recently, extended Lya emission was discov- 
ered around the z = 6.4 radio-quiet quasar CFHQSJ2329-0301 
(IWillott et alJ20lTtlGoto et "allBOLl . 

Extended Lya emissi on has also been observ e d around high- 
redshift radio galaxies (iHeckman et al.l 1 1991 at Ivan Oiik et al. | 
1996) , as well as around radio-quiet quasars (IChristensen et al.l 
2006: IBunker et alT2 003). The Lya flux from the nebula is about 



0.5% of the flux in the integrated broad Lya line of the QSO, in 
the case of radio-quiet quasars (RQQs). For radio-loud quasars 
(RLQs) or radio-gala xies, the Lyg flux of the nebula is an order 
of magnitude higher (Christ ensen et al.ll2006l) . presumably be- 
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cause the emission of RLQ gaseous envelopes is enhanced by in- 
teractions with the radio jets. In addition, it seems t hat the RLQs 
gener ally reside in richer environments than RQQs (lEllison et al.l 
ifool) . 

Two main power sources have been suggested to explain 
the luminosity of the Lya blobs (excluding the envelopes of 
RLQs). The first is photo-ionization b y a central source such 
as an AGN or a s tarburst region (e.g. IChelouche et ail 12008). 
iGeach et"ai1 J2009) recently presented arguments in favour of 
this idea, based on deep X-ray observations of Lyo- blobs. They 
unveiled obscured AGNs but found no diffuse X-ray emission 
that would have betrayed the existence of a hot gas component 
with T > 10 7 K; the lack of such a component rules out in- 
verse Compton scattering of CMB photons as an ionizing source. 
The other possible power source is cold accretion, a scenario 
where pristine gas falling into the potential well of a dark matter 
halo gets heated and ionized by collisions, co nverting gravita- 
tional energy into Lya radiation. According to Dijkstra & Loebl 
(2009), observed Lya blobs could be explained in this way if 
only > 10% of the gravitational binding energy of cold gas be- 
ing accreted onto a galaxy is conve rted into Lya radiatio n (such 
a conversion was also discussed by Haiman et al. 2000 and by 
Farda l et al. 1 120011) . In addition, since the emission would come 
from filamentary cold flows, it would not depend on the exis- 
tence of a central ionizing source such as an AGN, because the 
cold gas would be self-shielded from the ionizing radiation. 

A blind search for Lya blobs has the advantage of providing 
a bias-free sample, but is very time-consuming. It requires wide- 
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field, narrow-band imaging, and the candidates found in this 
way still have to be observed spectroscopically to be confirmed. 
Alternatively, looking for Lya blobs around radio-quiet quasars 
takes advantage of pre-existing catalogues. In addition, bright 
quasars are thought to reside in massive halos (~ 5 x 10 12 M Q , 
see e.g. lHopkins et a l. 2007), which are a lso thought to host the 
so-called "cold flows" (Keres et al. 2009) that could emit a sub- 
stantial Lya flux dDiikstra & Loebir2009l) . The main drawbacks 
of this alternative search technique are that (1) it is not guaran- 
teed a priori that a Lya nebula is present around each quasar, (2) 
the sample will inevitably be biased, and (3) the bright quasar 
emission has to be subtracted before one can see and character- 
ize the Lya- envelope. 

We initiated, a few years ago, a pilot survey aimed at ex- 
ploring in detail the spatial extent, the luminosity and kinemat- 
ics of the large hydrogen envelopes of remote quasars spanning 
a broad magnitude range, at redshift z ~ 4.5. In this way, the 
possible relation between the quasar luminosity and the extent 
and luminosity of the envelope could be explored. In a model 
where photoionizati on by the AGN is the main cause of the 
Lya luminosity (e.g. lHaiman & Reesll2001l) . one should expect 
the latter to be strongly correlated with the luminosity of the 
quasar, at least that coming fr om the ionizing Lyman c ontinuum. 
Alternatively, in the model of lDiikstra & LoebT d2009). no direct 
relation between the AGN luminosity and that of the blob is ex- 
pected. Preliminary results of our survey have been published 
dCourbin et al.ll2.Q08l Paper I) for three quasars, one of which is 
surrounded by a Lya envelope (it appears that the marginal de- 
tection of a second one was spurious); the observational strategy 
was explained in more details in that paper. The present arti- 
cle describes our detection of a Lya nebula around two more 
quasars, and a marginal detection around a sixth one. The obser- 
vations are described in Section|2]and the results are presented in 
Section[3] The results are discussed assuming a cosmology with 
H = 72 km s _1 Mpc" 1 , Q. M = 0.3, and Q. A = 0.7. 



Table 1. Journal of observations and main characteristics of the 
quasars. 



Name JD(start) 


Exposure 


Airmass 


Seeing 


-2400000 


time (s) 


(start) 


(") 


SDSSJ14472+0401, z 


=4.580, R(AB)=19.9±0.2, M R =-2S.2 


54647.617 


1300 


1.30 


0.54 


54647.633 


1300 


1.38 


0.57 


54653.565 


1300 


1.18 


0.61 


54653.581 


1300 


1.22 


0.60 


54654.562 


1300 


1.18 


0.85 


54654.577 


1300 


1.22 


0.88 


54654.600 


1300 


1.31 


0.89 


54654.615 


1300 


1.39 


0.90 


SDSSJ21474-0838, z= 


=4.588, R(AB)=18.7±0.2, M«=-29.3 


54645.875 


1300 


1.07 


1.04 


54645.891 


1300 


1.10 


0.99 


54647.700 


1300 


1.48 


0.93 


54647.715 


1300 


1.35 


0.90 


54647.736 


1300 


1.23 


0.69 


54647.752 


1300 


1.17 


0.91 


54647.772 


1300 


1.11 


0.55 


54647.787 


1300 


1.08 


0.57 


BR 2237-0607, z=4.550, R(AB)= 


18.3±0.2, M R = 


-29.8 


54653.750 


1300 


1.29 


0.69 


54653.766 


1300 


1.21 


0.68 


54653.790 


1300 


1.13 


1.02 


54654.716 


1300 


1.52 


0.97 


54654.732 


1300 


1.39 


1.14 


54654.751 


1300 


1.27 


1.26 


54654.767 


1300 


1.20 


1.18 


54656.803 


1300 


1.09 


0.82 


54656.818 


1300 


1.07 


0.91 



Notes. The apparent magnitudes are given in the AB system. They are 
computed by integrating the quasar's spectrum through the RSPECIAL 
ESO filter curve. The absolute magnitude assumes H = 72 km s" 1 
Mpc-'and (fl„, C1 A ) = (0.3, 0.7). 



2. Deep VLT optical spectroscopy 

2.1. Sample and observations 

The quasars were selected from the 5th edition of the SDSS cat- 
alogue of quasars dSchneider et al. 2007) an d from the 12th edi- 
tion of the Verons' catalogue dVeron-Cettv & Veronll2006l) . All 
six objects are listed in Table_QSO of the Verons' catalogue with 
an asterix before their name, which means that they have not 
been explicitly associated with a radio source in the literature. 
One of th em, BR1 033-032 7, was explicitly considered as radio- 
quiet bv lKellv et alj ( l2008t) . The other objects that were optically 
detected are probably radio quiet, though statistically a possibil- 
ity remains that one of them might be radio-loud. The fraction 
of radio-loud objects i n quasars emittin g close to the Eddington 
limit is indeed < 5% dRafter et al.ll2009l Fig. 9); thus, there is 
at most one chance in three that one of our six quasars is radio - 
loud. 

The observations of the three new targets were carried out 
in ESO Period 81 (April-September 2008) in service mode with 
the FORS2 multi-object spectrograph attached to VLT-UT1. The 
ESO grism G1200R+93 has a resolving power R = 1070 with 
a 2" -slit, which ensures that we can catch most of the flux of 
the Lya envelope. This grism is used in combination with the 
GG435+81 order separating filter, leading to the wavelength 
coverage 6000 A < A < 7200 A. The maximum efficiency of 
this combination coincides well with the expected wavelength 
of the redshifted Lya line, i.e. about 6686 A. 



The multi-slit MXU mode is used, with slits that are long 
enough (typical length: ~ 20") to reliably model and subtract 
the sky emission. Only one scientific target is observable in each 
field, but the MXU capability is used to observe several stars 
through identical slits. In this way, a spectral point spread func- 
tion (PSF) is measured simultaneously with the quasar. This is 
crucial for the spatial decon volution of the data to work effi- 
ciently (Courb in et al.ll20 00) and to clearly separate the quasar 
spectrum from that of the putative envelope. 

To properly remove cosmic rays, the total exposure time is 
split into several shorter exposures. The journal of the actual ob- 
servations is presented in Table [TJ 

2.2. Reduction and spatial deconvolution 

The data reduction was carried out using the standard IRAF pro- 
cedures. The individual spectra listed in Table[T]were flat-fielded 
using dome flats, and wavelength-calibrated in two dimensions 
in order to correct the sky emission lines for slit curvature . The 
scale of the reduced data is 0.76 A per pixel in the spectral di- 
rection and 0.25" in the spatial direction. 

The sky emission was then subtracted from the individual 
frames by fitting a second order polynomial along the spatial 
direction. This fit considered only ten pixels on each side of the 
slit, and the sky at the position of the quasar was interpolated 
using this fit, both on the quasar and the PSF stars. 
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possible to estimate the luminosity of the Lya emission "under- 
neath" the quasar. Subsampling of the data is also possible with 
the MCS algorithm, hence the pixel size in Fig. Q~|is half that of 
the original data, i.e., the new pixel size is 0.125". 



Iff- 



J21 474-0838 (z=4.588) 




Fig. 1. The three combined and sky-subtracted FORS2 spectra 
along with their spatial deconvolution. The spatial resolution af- 
ter deconvolution is 0.25". In each panel, the vertical arrow in- 
dicates the position of the Lya- emission line of the quasar, the 
wavelength of which is indicated. The horizontal double arrow 
is 50 A wide. A Lya envelope is detected in all three quasars. 
The envelope of SDSS J14472+0401 is blueshifted from that 
of the quasar, while that of SDSS J2 1474-0838 and that of 
BR 2237-0607 are redshifted. The height of the spectra is 16" 
in all panels. 



3. Results 

The spectra of the three newly found Lya envelopes are shown 
in Fig. [2] toget her with the spectra o f the corresponding quasars 
(see Fig. 4 of ICourbin et alJ (120081) for the three other quasars 
previously observed). The Lya luminosity of the envelopes, their 
angular size, and their mean surface brightness are presented 
in Tables [2] and [3] We integrated the deconvolved spectra in 
the FORS2 RSPECIAL filter, which is also used to obtain short 
acquisition images prior to the long spectroscopic exposures. 
These "spectroscopic" AB magnitudes are given in Table [T] We 
checked that they are compatible with the simple aperture pho- 
tometry obtained from the short images. 

A Lya envelope was detected in four objects out of six, 
which suggests that these envelopes are present in two thirds of 
the QSOs at z ~ 4.5. Our flux limit (integrated over the whole 
slit and Lya- spectral profile) is indicated in Table [2] and dis- 
cussed below. The lack of envelope around SDSS J0939+0039 
and Q 2 1 39-4324 may be real, but could also be due to an unfor- 
tunate orientation of the slit, if the latter is oriented perpendicular 
to the long axis of the envelope. 

The measurable extent of the envelopes varies (in terms of 
radius) from r ~ 20 kpc for BR 2237-0607, to r ~ 32 kpc 
for BR 1033-0327, measured from the quasar's centroid to the 
noise level in the spatial profiles shown in Fig. [3] 

The surface brightness of the Lya fuzz is unaffected by slit 
losses, but the total luminosity is. Assuming that the Lya en- 
velopes are uniform face-on disks with diameters equal to the 
extents quoted in Table [2] we can estimate the amount of flux 
missed by using a slit width of 2". The observed luminosities of 
the Lya envelopes for all four quasars are given in Table [3] as 
well as the luminosities after correction for the slit clipping. 



Th e cosmic rays were removed using the L . A . Cosmic algo- 
rithm (Ivan Dokkuml200ll) . All the frames were visually checked 
after this process to ensure that no signal was mistakenly re- 
moved from the data. 

The shape of the spectra along the spectral direction is 
slightly distorted, i.e., the position of the spectrum changes as 
a function of wavelength. These distortions are corrected for and 
the spectra are eventually weighted so that their flux is the same 
at a reference wavelength, before they are combined to form a 
deep two-dimensional (2D) spectrum. We show in Fig. Q] the 
combined spectra of the three new quasars, after binning in both 
the spatial and the spectral directions. An extended Lya enve- 
lope is already visible in all three objects. 

We spatial ly deconvolved the spectra following the method 
described in ICourbin et alJ (l2000l) . which is an adaptation to 
spectroscopy of the "MCS" image deconvolution algorithm 
(Mag ain et al.l ll998). The results of this deconvolution are dis- 
played in Fig.Q] The algorithm uses the spatial information con- 
tained in the spectrum of several PSF stars in order to sharpen the 
data in the spatial direction. At the same time, it also decomposes 
the data into a point-source and an extended-source channel. The 
output of the deconvolution procedure consists of two individual 
spectra, one for the quasar and one for its host galaxy (or Lya 
envelope), free of any mutual light contamination. It is therefore 



3. 1 . Detection limit 

We now discuss the detection limit. This discussion was not in- 
cluded in paper I, where the detection limit was too optimistic 
because it was based on the deconvolved images rather than on 
the original ones. We now redefine this limit for all six objects 
observed so far (three in Paper I, three in the present work, see 
Tables and E). We use the case of BR 2237-0607 for this dis- 
cussion, because it is the target with both the longest exposure 
time and the faintest envelope if any. 

Using spectrophotometric standards, we verified the reliabil- 
ity of the flux calibration by computing the /? spe ciai magnitude 
of the quasar from its observed spectrum, using the /? spe ciai filter 
transmission curve: the obtained magnitude proved to be consis- 
tent with that obtained from the imaging observations using the 
same filter. The maximum intensity of the quasar BR 2237-0607 
is about 370000 ADU (analogic digital units), for one pixel 
along the wavelength axis (representing 0.76 A), after summa- 
tion along the spatial axis. Because the corresponding physical 
flux is 8.41 x 10~ 16 ergs -1 cirT 2 A~' , this extracted pixel (i.e. 
the pixel of the extracted quasar spectrum) receives a flux of 
6.39xl0~ 16 ergs~' cirT 2 on its 0.76 A width. Conversely, 1 ADU 
corresponds to 1.73 x 10~ 21 ergs -1 cirT 2 in one pixel. We now 
consider an extended source rather than a point source: the 8.23 
pixels that cover the slit width (considering only one pixel in the 
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Fig. 2. Extracted ID spectra for the three quasars. In each case, the top panel shows the quasar alone and the bottom panel shows 
the Lya envelope alone, after spatial deconvolution of the spectra. The vertical dotted line indicates the position of the un-absorbed 
Lya emission line, at the redshift of the quasar (see subsections 3.1 and 2.1). In all cases, the spectrum of the Lye envelope has 
been smoothed using a boxcar of 8 A. 

Table 2. Main properties of the Lya envelopes for all six quasars measured to date (Paper I and this work). 



Object 



A FWHM Extent 

[A] [A] (";kpc) 



Integrated flux 3 - cr detection limit Surface brightness 
[erg s _I cirr 2 ] [erg s -1 enr 2 ] [erg s~' cirr 2 "~ 2 ] 



SDSS J0939+0039 










(7; 45) 






3.9 x 10" 


ita 


< 2.8 X lO -1 *" 


BR 1033-0327 


6725.0 


±0.5 


50 ± 


10 


10;64 


1.4(= 


=0.1) x 10~ 17 


6.7 x 10" 


18 


7.2(±0.5) x 10- 19 


SDSS J14472+0401 


6756.3 


+ 0.5 


30 ± 


3 


6;38 


4.3(= 


=0.7) x 10~ 18 


2.5 x 10" 


18 


3.6(±0.6) x 10- 19 


Q 2139-4324 


6641.0 


±0.3 






(7; 45) 






2.2 x 10- 


8a 


< 1.6 x 10- 19fl 


SDSS J21474-0838 


6808.2 


± 1.0 


50 ± 


8 


8;51 


3.1(= 


=0.2) x 10~ 16 


3.3 x 10" 


18 


1.9(±0.1) x 10~ 17 


BR 2237-0607 


6773.6 


±2.0 


19 ± 


2 


4;26 


2.6(= 


=0.8) x 10~ 18 


2.4 x 10- 


18 


3.3(±1.0)x 10- 19 



Notes. All parameters are given in the observer's frame. The 3 — cr detection limit is based on the standard deviation of the background noise 
according to equ. [2] The surface brightness is integrated in wavelength but given per arcsec 2 , while the 3 - cr limit is spatially and spectrally 
integrated. 

assuming an average extent of 7" 



spatial direction) correspond to a solid angle of 0.252" x 2" = 
0.504" 2 on the sky. Therefore, there are 8.23/0.504 = 16.33 pix- 
els within a square arcsecond. Taking into account an observed 
standard deviation of cr, ~ 40 ADU (for individual pixels) in 
the final 2D spectrum (corresponding to a total exposure time of 
1 1 700 s), one sees that the standard deviation in the total back- 



ground intensity, summed over a square arcsecond is (outside the 
sky emission lines) 

cr(l" 2 ) a VTr33xo- ; ~ 484 ADU ~ 8.4xlO~ I9 ergs-' cirr 2 .(1) 

This is also the error in a signal spread over 1 arcsec 2 , pro- 
vided that it is small compared to the sky background (which 
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Table 3. Lye luminosity of the quasars in the BLR, compared with the luminosity of the extended envelopes. 



Object L(BLR) L(Lya) LconectedCLya) 

[ergs -1 ] [ergs -1 ] [ergs -1 ] 



SDSS J0939+0039 


4.1(d 


=0.2) 


X 


10 44 


< 7 


4x 10 41 " 


< 2 


x 10 42fl 


BR 1033-0327 


7.2(= 


=0.4) 


X 


10 45 


2.7(4 


0.2) x 10 42 


(1.14 


0.1) x 10 43 


SDSS J14472+0401 


2.1(= 


=0.1) 


X 


10 45 


8.6(=t 


:1.4)X 10 41 


(2.0=t 


0.3) x 10 42 


Q 2139-4324 


4.6(= 


=0.2) 


X 


10 44 


< 4 


1 x 10 41n 


< 1 


1 x 10 42fl 


SDSS J2 1474-0838 


6.8(+0.3) 


X 


10 45 


6.2(4 


:0.4) X 10 43 


1.9(4 


0.1) x 10 44 


BR 2237-0607 


1.4(= 


=0.1) 


X 


10 46 


5.1(4 


.1.6) x 10 41 


8.0(4 


2.5) x 10 41 



Notes. The flux in the quasar BLR is measured in the wavelength interval 1 200- 1230 A (rest-fra me). The last column gives the Lya- flux of the 
envelope, after correction for slit-clipping (see text). The objects published in Courbin et al. (2008) are included for the sake of completeness. 
<o) assuming an average extent of 7" 



is the case here) and assuming that the latter is sufficiently well- 
defined to have a negligible error. Therefore, the 3-cr minimum 
signal that can be considered significant after summation across 
the whole extent of the nebula (spanning n arcsec 2 ) is 



S(3cr) 12.1 x V«xo-; =484 V« ADU. 



(2) 



As a second consistency check, we used the Exposure Time 
Calculator forFORS2 (ESO version 3.2.7) and were able to con- 
firm the above figures: the flux corresponding to 1 ADU is re- 
covered to within 30% (so the validity of our flux calibration 
appears to be robust), and the standard deviation in the synthetic 
2D spectrum, resulting from a single 11700 s exposure, is only 
slightly below that of the empirical one (~ 30 ADU instead of 
- 35-40 ADU). 

In view of the above estimates, it appears that in paper I, 
we were too optimistic about the detection limit in the case of 
Q 2139-4324: our present calculations imply that the envelope 
cannot be assumed to have been detected, and only an upper 
limit can be assigned to its flux. 

3.2. Comments on individual objects 

We now comment on briefly each of the three newly discovered 
objects. 



3.2.1. SDSS J1 4472+0401 

The extended emission is blue-shifted with respect to the ex- 
pected position of the Lya line of the quasar by AA = -20 + 4 A, 
i.e., AV = -889 + 180 kms -1 . The redshift of the quasar may 
be underestimated, as it is only measured from both the Lya line 
itself and the pair of CIV ,11549 and CIII] ,11909 lines, all of 
which are known to be blueshifted with re spect to the actual red- 
shift of the quasar (Mcintosh et al.lll999l) . The Lya line, on the 
other hand, is strongly affected by the Lya forest. If the red- 
shift of the quasar were underestimated, one would expect the 
surrounding Lya emission to be redshifted with respect to the 
expected position of the Lya line, while we see the reverse. In 
view of the strong spatial asymmetry of the emitting region (see 
Fig- E and [3}, it is possible that the "envelope" is not centred at 
all on the quasar, but is a blob at some distance away from it. 
The blue shift may be produced by radiative transfer eff ects in 
both a c ollapsing neutral gas clou d (Diikstr a et al.ll2QQ6allbb and 
shocks dNeufeld & McKedl 19881) . 

Interestingly, the nebular emission is matched by a deep ab- 
sorption feature in the quasar spectrum, at exactly the same 
wavelength. This is reminiscent of the proximate damped Lya 
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Fig. 3. Spatial profiles of the envelopes. The black, quasi-vertical 
lines show the intensity profile of the quasar, while the red line is 
the intensity profile of the envelope, integrated across the spec- 
tral width of the Lya line (i.e. an arbitrary range of 20 A for 
objects with no significant envelope, and a 1 x FWHM for the 
other objects). The blue and green lines show the noise, i.e. the 
profile integrated across wavelength ranges of the same width, 
but centred 76 A bluer, respectively redder than the envelope 
emission line (except for BR 2237-0607, for which the noise 
was sampled at ±29 A from the emission line). 



system (PDLA) discovered by iHennawiet aD ((2009) and asso- 
ciated with the quasar SDSS J124020.91 + 145535.6, except that 
the Lya absorption does not seem to be damped in our case. 
Some non-negligible flux remains at the bottom of the line, and 
there is no indication of Lorentzian wings in the absorption pro- 
file. Nevertheless, the width of this absorption is large, with 
FWHM = 2100 + 100 kms -1 , which is about 50% wider than 
the emission line of the blob; the lack of damped wings suggests 
that the FWHM corresponds to the velocity dispersion of the ab- 
sorbing gas. On the other hand, the emission line seems to have a 
double structure, which, if taken into account, reaches about the 
same width as that of the absorption line. Higher dispersion and 
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higher signal-to-noise ratio spectra of this object would allow us 
to clarify the latter point. 

3.2.2. SDSSJ21 474-0838 

This is by far the quasar with the brightest Lya blob in our sam- 
ple. It is an order of magnitude brighter than the second bright- 
est one, which surrounds the quasar BR 1033-0327. As in the 
latter object, the emission is slightly redward of the expected 
position of the Lya line of the quasar, by AA = +32 + 4 A or 
AV = +1400 ±180 km s -1 . This is very close to, though slightly 
larger than, the shift measured in BR 1033-0327, and the expla- 
nation of it might be the same, i.e. the known blueshift of the 
Lya, C IV A1549, and C III] /11909 lines. Interestingly, the emis- 
sion line appears to be slightly asymmetric, with a sharper drop 
on its red side than the rise on its blue side; a similar asymmetry 
is seen in the emission line of BR 1033-0327. This asymmetry 
might also be present in SDSS J 14472 +0401, though it is less 
prominent and possibly due to an absorption feature on the blue 
side of the emission line. 

The positive flux seen on Fig. [2] on both sides of the enve- 
lope emission line (and especially on its red side), is probably 
spurious. There is a bad column on the CCD, running parallel to 
the quasar spectrum, about 28 pixels away from it (or 7"), which 
may have caused problems with the sky subtraction. The quasar 
is so bright that a small imperfection in the PSF profile might 
have produced a pseudo-continuum. 

3.2.3. BR 2237-0607 

This quasar shows the faintest detected Lya- blob in our sample 
(disregarding the two non-detections), with a luminosity about 
two times lower than the second faintest blob in our sample 
(SDSS J 14472+0401). This is also the most difficult case, since 
the quasar is the brightest in the sample, so that the contrast be- 
tween the quasar and the envelope is the largest. This is why 
it is detected only at the ~ 3 - <x level, according to the crite- 
rion presented in Section 3.1. The emission is strongly redshifted 
relative to the expected position of the Lya- line of the quasar: 
AA = +87 ± 5 A or A V = +3850 ± 220 km s" 1 . Such a consider- 
able shift cannot be explained by absorption, since the Lya line 
of the quasar does show significant flux down to at least 6730 A 
while the nebula emits only between ~ 6760 A and ~ 6790 A. 
Thus, one has to admit that either the redshift has been underesti- 
mated, or the nebula is receding at high velocity from the quasar; 
it is unlikely that radiative transfer effects could mimick such a 
fast motion, since one would require a high absorption/emission 
cross-section in the line profile at these velocities. 

3.3. Redshift dependence 

The observed surface brightness of the four Lya- envelopes de- 
tected in our programme is, on average, fainter by about 1-2 
orders of magnitude than for the CJW objects. Our brightest 
envelope has about the same surface brightness and integrated 
flux as typical CJW envelopes, while our faintest ones are a 
hundred times dimmer. This contrast cannot be ascribed exclu- 
sively to the different redshift ranges of the two samples. Our 
sample is at z\ ~ 4.5, while most of the objects in CJW are at 
Zi ~ 3.3, so the flux ratio by redshift dimming alone would be 
F (Z2)/F(zi) = (1 + z 2 ) 4 /(l + zi) 4 ~ 2.7, i.e., much smaller than 
the observed ratio. 



Therefore, taken at face value, our sample as a whole sug- 
gests that higher-z objects tend to have lower surface bright- 
nesses. The total luminosities of the envelopes in our sample 
overlap those in CJW: SDSS J21474-0838 has a higher Lya 
luminosity (when corrected for slit clipping) than all objects in 
CJW, the envelope of BR 1033-0327 has a luminosity close to 
the average one in the CJW sample, and the two other objects are 
fainter. Thus, the envelope luminosities in our sample seem to 
more closely match those of CJW than the surface brightnesses 
do, but this may be due to our deeper sensitivity. Observations 
of the CJW envelopes with the same sensitivity as ours would 
probably reveal that they have faint extensions. Consequently, 
the envelopes would be found to have higher total estimated lu- 
minosities, on the one hand, and a lower average surface bright- 
ness, on the other. The spatial extent of the envelopes in our sam- 
ple ranges from ~ 26 kpc to ~ 64 kpc, while that of CJW ranges 
from ~ 10 kpc to ~ 60 kpc. That explains why the luminosities 
of our envelopes are often similar to those of CJW, in spite of a 
much lower surface brightness. The mean size of the envelopes 
in our sample is indeed r ~ 45 kpc compared to r ~ 26.4 kpc in 
CJW's sample. This translates into a ratio of 2.9 in area. 

The quasar luminosities are slightly larger in CJW's sam- 
pl e than in ours. The ave r age absolute magnitude Mb, as listed 
in IVeron-Cettv & Veronl (120061) . is -29.0 for CJW's complete 
sample of seven quasars, while it is -28.2 for our sample of six 
quasars. Since a correlation seems to exist between the luminos- 
ity of the quasar and that of the Lya envelope (see below), such 
a magnitude difference may explain part of the difference in Lya 
envelope luminosities. 

A sample of quasars should be observed at z ~ 3.3 with the 
same technique and depth as used at z ~ 4.5 in order to settle the 
question. 

3.4. Dependence on the luminosity of the quasar 

The total flux in the Lya envelope is shown in Fig. [4] as a func- 
tion of the quasar luminosity in the broad Lya line. After cor- 
recting for slit-clipping, we find that most points, combined with 
the ones of CJW, follow the linear relation 



log 



LiLyo) 



UBLR) 



= (0.65 + 0.53) x log[L(BLR)] - (32.3 ± 24.5) ,(3) 



suggesting that L(Lya) varies strongly with L(BLR). The regres- 
sion line was determined after excluding BR 2237-0607 (the 
two objects with upper limits to L(Lya) were not considered). 
The scatter is 0.39 dex and the correlation coefficient is 0.874, 
leaving no doubt about the statistical significance. Nevertheless, 
BR 2237-0607 remains two orders of magnitude below the lin- 
ear relation, which seems difficult to explain by slit-clipping ef- 
fects alone. Since SDSS J0939+0039 and Q 2139-4324 may 
also lie well below the regression line, the latter should perhaps 
be considered as an upper envelope rather than as a real one- 
to-one relation proper. However, the larger scatter of our points 
compared to those of CJW around the regression line, may sim- 
ply reflect our larger uncertainty in L(Lya) due to the slit clip- 
ping and unknown shape of the envelope. CJW did not have this 
problem because they used integral field spectroscopy. 

The linear fit implicitly assumes no strong redshift evolution 
of the luminosity of the envelopes, since it relies on a mix of 
objects at redshifts z ~ 3 and z ~ 4.5. Our new observations for 
the three objects therefore seem to support the trend that brighter 
quasars also have brighter Lya envelopes, under the assumption 
of negligible redshift evolution. 
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Fig. 4. Top: Relation between the total luminosity of the Lya en- 
velope and that of the quasar in the broad Lya line. Our measure- 
ments (red squares) are compared with CJW (black circles). The 
open squares represent the direct measurements, while the filled 
dots are those corrected for slit clipping (see text). The down- 
ward arrows show our upper limits for SDSS J0939+0039and 
Q 2139-4324, where no Lya envelope is seen. In computing 
this limit, we assumed a radius r=23 kpc (i.e., the mean size of 
the other four objects) and we did not correct for slit-clipping. 
The black solid line was least squares fitted to the objects of 
CJW and to our detected envelopes (with correction for slit clip- 
ping), except for BR 2237-0607. The blue triangles relate the 
Lya luminosity (corrected for slit clipping) with the luminosity 
of the ionizing radiation from the UV continuum of the quasar 
(A < 912 A, see Subsection 3.3). The arrows again show the up- 
per limits to the luminosity of the undetected envelopes. Bottom: 
The radius of the envelope as a function of L(BLR). Note the 
tight trend followed by the envelopes of the bright RQQs of 
CJW. Note also that BR 1033-0327 stands above the maximum 
size of the RQQs of CJW, likely owing to our deeper flux limit. 



The relation between the size of the envelope and the quasar 
luminosity is much less clear. CJW find that brighter envelopes 
tend to be larger, which, combined with the above correlation 
between envelope and BLR luminosities, implies that brighter 
quasars should be embedded in larger envelopes. We find no 
such trend on the basis of our results for our four objects, as 
shown in Fig. [4] (bottom panel). The marginal trend obtained 
with the CJW data may be explained by the difficulty in sub- 
tracting the quasar light: the brighter the quasar, the wider the 
envelope has to be in order to be detected unambiguously. Our 
deeper observations, coupled with a cleaner subtraction of the 
quasar light, make our results less sensitive to this effect. In ad- 
dition, the small field of view used in CJW implies that there has 
been some severe clipping of the envelopes, if they extend much 
beyond a few arcsecs. 

We now examine the possible correlation between the ion- 
iz ing flux of the quasar a nd the Lya flux of the envelope. As 
in Christe nsen et al.l d2006l) . we use the template spectrum based 
on the Sloan Digital Sky Survey (IVanden Berk et al.l20 01). com- 
pleted by th e composite Far U ltraviolet Spectroscopic Explorer 
spectrum of IScott et al.l d2004l) . For simpli city, we adopted the 
power-law fit F v <x y-°- 56 determined by IScott et al] (120041) to 
represent the quasar spectrum in the far ultra-violet, from 800 A 
to 912 A. We redshifted this composite template spectrum to 
the observed value of each of our observed targets, and adjusted 
its intensity so that it matched the observed flux in the range 
~ 1286 - 1291 A in the rest frame (~ 7071 - 7099 A at z = 4.5). 
Then, we integrated the spectrum between 800 A and 912 A af- 
ter having rebinned it to the rest frame, obtaining the ionizing 
luminosity L(< 912A). This luminosity depends of course on 
the lower integration limit, which is arbitrary. Another choice 
would have simply changed L(< 912A) by a constant factor. 
Figure [4] (upper panel) shows the Lya luminosity of the enve- 
lope (corrected for slit clipping) as a function of the ionizing 
luminosity defined above. At first sight, there is no correlation, 
which tends to confirm the result obtained by Christ ensen et al.l 
(2006), but at larger redshifts (z ~ 4.5) and for a wider range 
of Lya luminosities: about 2.3 dex instead of ~ 1 dex, taking 
only the RQQs into account. However, the lack of correlation 
is essentially due to BR 2237-0607, and removing this object 
would result in a clear trend, similar to the one found using the 
BLR luminosity instead of L(< 912A). Hence, the interpretation 
of our results depends on the weight given to BR 2237-0607: if 
one fully takes it into account, then there is no correlation and 
this could be interpreted as an argument against the quasar being 
the main cause of the Lya emission of the envelope; conversely, 
excluding BR 2237-0607 as a pathological case restores the cor- 
relation and leads to the opposite conclusion. In any case, small 
number statistics and slit-clipping effects make such a conclu- 
sion fragile, as well as the uncertainty inherent to the use of a 
uniform template spectrum for the quasar. 

3.5. Surface brightness and width of the emission line 

The FWHM of the Lye emission line varies by only a factor of 
~ 2.5 from one object to the other in our sample, but it may 
be interesting to consider whether it is correlated with e.g. the 
Lyor luminosity of the nebula. The latter spans a wide range, of 
about 2.3 dex. Plotting log(FWHM) against log L(Lya), we in- 
deed found a rough trend of increasing FWHM (corrected for the 
intrumental width) with increasing luminosity, although there 
are four points only, so it is not statistically significant. The 
correlation coefficient of 0.76 and Student t test of 2.03 leave 



North et al.: Extended Lya envelopes around z = 4.5 quasars 



about a 20% probability of finding the observed correlation by 
chance alone. Furthermore, the CJW data do not fit the correla- 
tion, which confirms that it cannot be real. 

Figure [5] presents the average surface brightness of the en- 
velope versus the average FWHM of its Lya line, for our four 
objects as well as for the RQ Qs studied by CJW , for th e quies- 
cent halos o f radio galaxies oflVillar-Martfn et aI1(l2003l) . and for 
the RLQs of iHeckman et al.l (Il991alfbl) . The surface brightnesses 
taken from the literature were all corrected to a common redshift 
of z — 4.54 (the average value of our sample), so that they may 
be compared. The FWHM values listed in Table 1 of CJW are 
local values, which do not include the velocity field of the gas; 
in order to compare with our values, which include a system- 
atic velocity field if any, we estimated "global" FWHM values 
from the one-dimensional (ID) spectra shown in Fig. 1 of CJW, 
which are larger since they inclu de the velocity field. Th e sur- 
face brightnesses of the objects of Heckman et ajj (Il991al) were 
estimated from their Table 1 , by dividing the spectral flux F s by 
1.5 x D s (1.5" being the width of the slit and D s the maximum 
angular extent of the envelope) when the envelope is detected on 
one side of the quasar only, and by dividing F s by 3 x D s when 
the envelope is detected on both sides of the quasar. 

It seems that the envelopes of all quasars (whether radio 
loud or not) and radio galaxies have similar maximum surface 
brightnesses. Most of our objects lie below the general trend, 
thanks to our much deeper detection limit, and because not all 
RQQs are surrounded by bright blobs. On the other hand, the 
FWHM of the extended Lye emission spans almost a factor of 
ten. Surprisingly, the object in our sample that has the highest 
surface brightness, has a wide FWHM, much wider than those 
of CJW, and surrounds a RQQ. Our objects with fainter surface 
brightnesses have, on average, wider emission lines than those 
of CJW. The reason for this is unclear and this result may be 
a statistical fluke. Nevertheless, it is interesting to note that the 
lower left region of the diagram is observationally more accessi- 
ble than the lower right region, i.e. the detection of low surface 
brightness envelopes is easier for small FWHMs than for large 
ones. Therefore, the relative scarcity of objects in this region 
cannot be the result of an observational bias, and must be real. 
One can only speculate that the kinematics of these envelopes 
might be more violent at redshift 4.5 than later, though such a 
rapid evolution (recalling that the CJW objects are at z ~ 3.2) 
seems doubtful. 

3.6. Upper limit to the flux in the A/v doublet 

The N V/H238.81 - 1242.80 A doublet is not detected in emis- 
sion, in any of the five envelopes detected here in Lya. It may 
therefore be interesting to estimate the upper limit to the flux 
emitted by this doublet, because this may be translated into an 
upper limit to the metallicity. We performed this estimate in the 
following very simple way. We adjusted the flux so that the dou- 
blet is almost detectable in the ID spectra of Fig. [2] under the as- 
sumptions that: 1) the width of the Lya line is essentially due to 
a Maxwellian distribution of the gas velocities, so that the same 
width can be adopted for each component of the doublet, 2) the 
two components of the doublet have the same strength, and 3) the 
intensities of the two components can simply be added. We are 
aware that our three assumptions are simplistic, especially the 
second one, which is equivalent to assuming an optically thick 
gas. The weighted oscillator strength gf of the blue component 
of the doublet is t wice as large as that of the red component 
dMarsh et al.lll995l) . so that the ratio of the blue to red compo- 
nent intensities would be two in the optically thin case. The lat- 
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Fig. 5. Average surface brightness of the envelope, as a function 
of the width of the Lya emission line (expressed in km s _1 ). The 
red dots represent our results. The surface of the points is pro- 
portional to the size of the envelope. 



Table 4. Upper limits to the flux emitted by the envelope in the 
N v doublet, compared to the flux emitted in the Lya line. 



Object 


Flux(Nv) 
[erg cm" 


Flux(Lyo') 


Flux 
ratio 


BR 1033-0327 


< 2x 10- 1S 


1.4 x 10-" 


< 0.14 


SDSS J14472+0401 


< 1 x 10- 18 


4.3 x 10~ 18 


< 0.23 


SDSS J21474-0838 


< 5 x 10~ 17 


3.1 x lO" 16 


< 0.16 


BR 2237-0607 


< 1.3 x 10~ 18 


2.6 x 10~ 18 


<0.50 



Notes. The last column gives the ratio of column 2 to column 3 (see 
text). 



ter case would make the detection easier because the resulting 
peak would be narrower; conversely, our assumption implies a 
wider peak, hence a less easy detection, and so appears rather 
conservative. 

The results are summarized in Table |4] which lists respec- 
tively the upper limit to the flux possibly emitted in the N v dou- 
blet, the flux in the Lye line, and the ratio of the two. The ratios 
are not very compelling, since the lowest is 0.14, meaning that 
the N v doublet is at least seven times less strong than the Lya 
line. One sees that, the narrower the Lya line, the less compelling 
the upper limit on the N v doublet, because the flux of the latter 
is distributed over a wider wavelength range. 

4. Conclusions 

We have performed deep slit spectroscopy of six radio-quiet 
quasars in the redshift interval 4.460 < z < 4.588 and found 
extended Lya emission around four of them. The depth of our 
detection limit is unprecedented, so we have been able to detect 
nebulae that are much fainter than in previous surveys. 

The main conclusions of this study can be summarized as 
follows: 
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- At z ~ 4.5, extended Lya envelopes are found around 
roughly two-thirds of the quasars. 

- The size of the typical envelope is very large (between ~ 
26 kpc and ~ 64 kpc), and does not depend on the Lya lumi- 
nosity of the BLR of the quasar. 

- The average surface brightness of the envelopes is 
very low, since it ranges from ~ 3 x 10~ 19 to 2 x 
10~ 17 erg s -1 cirT 2 arcsec -2 . Such low values are due to the 
large size of the envelopes. Both the large size and the low 
surface brightness seem difficult to reconcile with a scenario 
that attributes most o f the Lya luminosity to fluorescence in- 
duced by the quasar (Hai man & ReesluOOll) . 

- The luminosity of the envelopes seems to be correlated with 
that of the BLR, but this is not an absolute rule, because 
of the exceptional behaviour of BR 2237-0607. Likewise, 
the envelope luminosity appears to be roughly correlated 
with the predicted ionizing flux of the quasar, but only when 
BR 2237-0607 is excluded. If confirmed, this trend would 
resemble that of CJW for RQQs, but over a wider range of 
luminosities. 

- The average FWHM of the extended emission line varies 
from about 900 kms~' to 2500 kms -1 in our sample 
alone, while it is often smaller in other samples (do wn to 
250 km s~'in the radio galaxies o flwii-Martm 2007|). 

- The N v/11238.81 - 1242.80 A doublet remains undetected 
in all our objects. The most compelling non-detection occurs 
in BR 1033-0327, where the flux in the doublet is at least 
seven times lower than in the Lya line. 

D o these results favour the "cold accretion" scenario advocated 
bv lDijkstra & Loeblfc oog). or rat her the "heating" scenario sup- 
ported by e.g. Geac h et al.l d2009)'iPl As discussed above, our re- 
sults remain ambiguous, because of the unexpected behaviour 
of BR 2237-0607. Arbitrarily discarding this object, the pos- 
itive trend between the ionizing flux and envelope luminos- 
ity might be considered as favouring the heating scenario, but 
taken at face value our results do not provide enough evidence 
for either scenario. A larger sample would enable us to settle 
the issue, because a clear correlation between the AGN UV 
luminosity and the Lya luminosity of the envelope is not ex- 
pected in the cold flow model. In any case, future theoreti- 
cal works will have to t ake the observed low s urface bright- 
nesses into account. As Haiman & Reesl (1200 ll) stress in the 
last sentence of their paper: "While a detection of Lya fuzz 
would provide a direct probe of galaxy formation, nondetec- 
tions at the level of 10~ 17 ergs -1 cm -2 arcsec -2 would already 
imply strong constraints." Here, we have provided two non- 
detections, but at a deeper level (assuming the envelope does 
extend over at least a few arcsec 2 ), and three detections of 
envelopes with surface brightnesses significantly lower than 
10~ 17 erg s _I cirT 2 arcsec -2 . Only one object was found to have 
a surface brightness of the order of 10~ 17 erg s _1 cm -2 arcsec -2 . 

A third possibility might be that the detected emission is not 
actual emission but rather scattering of broad emission-line pho- 
tons from the quasar. This is suggested by the observed rela- 
tion L(Lya) oc L(BLR) L65±(l5i (Eq.rj, which is compatible to 
L(Lya) oc L(BLR) 2 . The scattered flux is, to zeroth order, pro- 

1 The "heating" versus "cooling" alternative highlighted by the title 
of Geach et al.'s paper is confusing. Indeed, in both models the cold 
gas is heated, which makes it able to glow in Lya. Only the source of 
heating differs. In one case, heating is provided by the UV radiation of 
a quasar, while in the other case heating results from the conversion of 
gravitational energy into thermal energy; in the first case, there is photo- 
ionization, while in the second case there is collisional excitation. 



portional to the BLR luminosity times the amount of matter in 
the envelope; if the latter scales with the bulge mass, which is 
itself proportional to the supermassive black hole mass, which 
determines the quasar luminosity, then it is also proportional to 
the BLR luminosity, indeed implying that L(Ly<?) oc L(BLR) 2 . 
Photons are then scattered by cool gas on large scales before 
reaching the observer. Despite the velocity shift between the gas 
and the quasar redshift (of typically a few hundreds of kms -1 ), 
there is enough photon flux from the broad quasar emission line 
(whose width is several thousands of km s -1 ) to produce a nebu- 
lar fuzz. Admittedly, this argument may prove excessively naive, 
e.g. because in the cold flow model, the amount of cold g as in- 
creases less rapidly than the halo mass dKeres et al.ll2009l) : the 
predicted power-law index depends on largely uncertain prop- 
erties of quasar environments at z — 4.5, but a more detailed 
treatment of these issues is beyond the scope of this paper. As 
for timescales, the nebular size is < 1 x 10 6 light years, hence 
much shorter than the typical quasar lifetime, which is thought 
to be of the order of 1 x 10 7 years. That the brighter the broad 
emission-line emission from the QSO, the stronger the nebular 
emission, is a result that is consistent with this scattering sce- 
nario. The material responsible for the scattering may originate 
in tidal tails and other debris, associated with the presumably 
interacting system leading to the onset of quasar activity. This 
would explain the asymmetry (of the order of 1 -2", correspond- 
ing to s ome ten kpc, see Fig. [3]l and the extent o f the emission 
(see e.g. lMihos & Hernquistll9 96; Springel 2000). The observed 
velocity difference between the peak emission of the QSO and 
the nebula, < 1000 km s -1 , is rather extreme for a local merger. 
Nevertheless, the quasar environments at z ~ 4.5 may be more 
extreme than those found at lower redshifts. These quasars may 
live in proto-cluster like environments with velocity dispersions 
of order 1000 kms -1 , which resemble today's richest clusters. 
The presence of an additional starburst wind component could 
contribute to the signal, as well as the cold gas being accreted 
onto the host galaxy (i.e. the cold flow gas). 

Additional spectroscopic measurements at the same redshift 
would be welcome in order to confirm the correlations found 
here. However, the crucial data that is badly needed to make 
progress in understanding these envelopes is narrow-band imag- 
ing, which will allow us to obtain much more reliable envelope 
shapes, extents, and luminosities. The latter are crucial discrimi- 
nating criteria to identify the mechanism responsible for the en- 
velope emission. 
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